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1. Introduction 



I ' universe. We compare the temperature correlation function and the angular power 

spectrum for the Poincare dodecahedral universe with the first-year WMAP data 
' and find that this multi-connected universe can explain the surprisingly low CMB 

anisotropy on large scales found by WMAP provided that the total energy density 
parameter f2tot is in the range 1.016 . . . 1.020. The ensemble average over the primordial 
' perturbations is assumed to be the scale-invariant Harrison- Zel'dovich spectrum. The 

■ circles-in-the-sky signature is studied and it is found that the signal of the six pairs of 

[ matched circles could be missed by current analyses of CMB sky maps. 

o 
o 

^ ■ The question of tlie global geometry of the Universe, i. e. of its spatial curvature, 
topology and thus its shape, is a very fundamental one. While the problem of the 
spatial curvature is discussed in all textbooks on cosmology, mainly in connection with 
' the total mass/energy densities in the Universe and the various inflationary scenarios, 
the possibility of a multi-connected Universe is only rarely mentioned. This is surprising 
in view of the fact that it was already realized before Einstein's seminal first paper on 
cosmology that the three-dimensional space of astronomy might be not only non- 
Euclidean, but also multi-connected according to a given Clifford-Klein space form (see 
e.g. |2]). If the Universe is assumed to be simply connected, the cosmological principle 
implies that it is not only locally, but also globally isotropic and homogeneous. This 
means that all spatial points are geometrically equivalent. In the case of a multi- 
connected Universe, however, the cosmological principle holds in general only on the 
universal covering space of constant curvature, i. e. on and for a fiat, positively 

and negatively curved Universe, respectively. While the Einstein gravitational field 
equations still hold in a mult i- connected Universe, the global structure of the Universe 
at large scales will be more complicated. 

The discovery of the temperature fluctuations 6T of the cosmic microwave 
background radiation (CMB) by COBE in 1992 ,3, and the detailed measurements 
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by WMAP j3] and by other groups have led to a renewed interest in the question of the 
global geometry of the Universe (see OE] for reviews). 

The most promising signatures to detect a possible multi-connected spatial 
structure of the Universe are the strange suppression of the CMB quadrupole and 
octopole, and of the temperature correlation function first observed by COBE j2] and 
nicely confirmed by WMAP j4j, and the so-called circles- in-the-sky-signature proposed 
in 0. 

Until about 1999, all observational data were consistent with the fact that the 
total energy density Stot (at the present epoch) amounted only to ca. 30% of the 
critical energy density ^crit = ^g^Q , i-e. fitot := £^tot/£^crit — 0.3 < 1, and thus 
strongly indicated that the geometry of the Universe is hyperbolic. Consequently, several 
groups performed detailed studies of compact and non-compact models of the Universe 
possessing hyperbolic geometry. Although later observations on the magnitude redshift 
relation of the supernovae of type la CD] provided evidence for a large dark energy 
component, and the determination of the first acoustic peak in the CMB anisotropy 
by TOCO [m Hg, BOOMERanG ^ and MAXIMA-1 O CSl indicated a nearly 
fiat Universe, i.e. f2tot — 1, the available data are still compatible with the spatial 
geometry of the Universe being hyperbolic. (For a discussion of hyperbolic universes, 

see [iniii7i[iHi[inii2niinii2ai23ii2iii2ai2nii2ii) 

The concordance model of cosmology assumes an exactly (spatially) flat Universe 
with the topology of and, furthermore, that the dark energy is given by a positive 

COSmological constant A, i. e. I^A := ^ = 1 - ^mat - ^rad with f^mat = ^bar + ^cdm, 

where the various parameters denote the present values of the baryonic (bar), cold 
dark matter (cdm), matter (mat) and radiation (rad) energy densities in units of ^crit 
(ACDM model). In the following, we shall refer to three variants of the concordance 
model presented by the WMAP team on the Legacy Archive for Microwave Background 
Data Analysis (LAMBDA) Web site http://lambda.gsfc.nasa.gov , see also |2H|. These 
models give a good overall fit to the CMB anisotropy on small and medium scales, but 
there remains a strange discrepancy at large scales as first observed by COBE [7J and 
later substantiated by WMAP 

Let us expand the temperature fluctuations 6T{n) of the microwave sky, where the 
dipole contribution has been subtracted, into real spherical harmonics Yim{n) on 5^, 

oo I 

5T{n) := X] aimYimin) , (1) 

1=2 m=-l 

where n denotes the unit vector in the direction from which the photons arrive. (Note 
that the monopole and dipole terms, / = 0, 1, are not included in the sum (0).) From 
the real expansion coefficients aim one forms the multipole moments 

m=—l 
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and the angular power spectrum 

--'^c. . (3) 

The average (...) in Q denotes an ensemble average over the primordial perturbations 
to be discussed in section 01 respectively an ensemble average over the universal 
observers. The temperature two-point correlation function C{'d) is defined as C{'d) : = 
{6T{n)6T{n')) with n- h' = cos^, which can be computed from the multipole moments 
(121) under the assumption of statistical isotropy as 

C(^9) ^ — 5^(2/ + l)GPKcos^) . (4) 

1=2 
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Figure 1. C(i?) from first year WMAP data (solid curve) in comparison to the three 
concordance models provided by the WMAP team at http://lambda.gsfc.nasa.govj 
model pi assumes a power law fit to the WMAP, CBI and ACBAR data, model ri 
assumes a running index fit to the WMAP, CBI and ACBAR data, and finally, model 
bf assumes a running index fit to the WMAP, CBI and ACBAR data and additionally 
takes the large scale data fi:om 2dF and Ly-alpha into account. 

In figure Ewe show C{'d) as solid curve as measured by WMAP |4j. The important 
conclusion to be drawn from figure^is that the temperature correlation function displays 
very weak correlations at wide angles, 60° ^ ^ ^ 160°. The data are compared with 
three different versions of the concordance model obtained by the WMAP team as best 
fits of their data combined with various other data. The model pi (dotted curve) assumes 
a power law fit, i. e. with no running index, to the WMAP j^, CBI [231201 and ACBAR 
pi] data, the model ri (dashed curve) assumes a running index fit to the WMAP, CBI 
and ACBAR data, and finally, the model bf (dashed-dotted curve) assumes again a 
running index fit to the WMAP, CBI and ACBAR data and additionally takes the 
large scale structure data from 2dF [22] and Lyman-a [HH] into account. It is seen from 
figure n that the difference between the three concordance models is not considerable. 
One makes, however, the important observation that the three concordance models 




Poincare Dodecahedron 



4 



using the best-fit values for tlie cosmological parameters as obtained by WMAP do not 
reproduce tlie experimentally observed suppression of power at wide angles, leaving 
us with "the mystery of the missing fluctuations" In the case of the largest 

angles, 160° < "(9 < 180°, where the WMAP curve in figure ^ shows even negative 
temperature correlations, one might argue that the measurements in the nearly back- 
to-back configuration are very subtle due to the need to apply a Galactic cut to the 
data and, furthermore, due to geometrical reasons, since the expectation values of C{i)) 
are computed from a much smaller number of pixels, and thus the discrepancy at these 
angles should not be overemphasized. There is, however, also a discrepancy between 
the three concordance models and the WMAP data at small angles, < 40°, which is 
clearly visible in figure Q 

Depending on certain priors, the WMAP team reported jlHj for the total energy 
density Qtot = 1-02 d= 0.02 together with fibar = 0.044 ± 0.004, Q^a^t = 0.27 ± 0.04, 
and h = 0.71^^^^ for the present day reduced Hubble constant (the errors give the la- 
deviation uncertainties). Taken at their face value, these parameters hint to a positively 
curved Universe. Recently Luminet et al. 35^ studied the Poincare dodecahedral space 
(for details, see section ^ which is one of the well-known space forms with constant 
positive curvature (see also [33 EZl Ell)- The authors of ref. computed the CMB 
multipoles for / = 2, 3 and 4, fitted the overall normalization factor to match the WMAP 
data at / = 4, and then examined the prediction for / = 2 and Z = 3 as a function of 
fltot- For ^tot = 1-013 they found a strong suppression of the power at Z = 2 and weak 
suppression at / = 3 in agreement with the WMAP data. Since the eigenfunctions of the 
Poincare dodecahedron are not known analytically, they must be computed numerically 
[38| l39] . In ref. jHS|, only the first three modes with wave number j3 = 13,21 and 
25 (comprising in total 59 eigenfunctions) have been used, which in turn restricted 
the discussion to the multipoles I < 4. There thus remains the question about how 
this extremely low wave number cut-off affects the prediction of the multipoles, since 
experience shows that increasing the cut-off usually enhances the integrated Sachs- Wolfe 
contribution. 

In this paper, we present a thorough discussion of the CMB anisotropy for the 
dodecahedral space topology which is based on the computation of the first 10521 
eigenfunctions corresponding to the large wave number cut-off f3 = 155. Taking within 
the tight-coupling approximation not only the ordinary, but also the integrated Sachs- 
Wolfe and also the Doppler effect into account, we are able to predict sufficiently 
many multipole moments such that a detailed comparison of the dodecahedral model 
with the WMAP data can be performed. We find that the temperature correlation 
function for the dodecahedral universe possesses very weak correlations at large scales 
in nice agreement with the WMAP data for fitot in the range 1.016. . .1.020. Besides 
the suppression of large-scale fluctuations, the dodecahedral model predicts 6 pairs of 
matching circles in the sky possessing an angular radius of 40° . . . 50° for the slight 
curvature of space corresponding to the cited range of fltot- 

We begin in section |21 by describing the spherical space forms and, in particular, the 
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Poincare dodecahedron, and discuss the general properties of the vibrational modes. In 
section El we outline our numerical method to determine the eigenmodes numerically. 
In section |31 we calculate the angular power spectrum and the temperature correlation 
function for the dodecahedral space and compare the predictions with the WMAP data. 
As a quantitative test, we study the S statistic introduced in For this statistic, it is 
found that nice agreement with the WMAP data is obtained for fitot = 1-016 . . . 1.020. 
In section we study the circles-in-the-sky signature using a quantity S introduced in 
|Hj (called S" in |HI). In the last section, there is a discussion and summary. 

2. The spherical space forms and their vibrational modes 

The three-dimensional spaces Ai^ of constant positive curvature K = +1 (spherical 
spaces) were classified by 1932 [101 HI] and are given by the quotient = S^/T of the 
3-sphere under the action of a discrete fixed-point free subgroup F C SO (4) of the 
isometrics of S^. All these manifolds are compact and are, apart from their universal 
covering S^, multi-connected. 

The unit 3-sphere is defined as the three-dimensional hypersurface of the 
four- dimensional unit ball in flat four- dimensional space with Cartesian coordinates 
{w, X, y, z) G M^, i. e. is given by the set of points which satisfy the condition 
w"^ + x"^ + y"^ + z"^ = 1. The spatial distance \dx\ between neighbouring points is 
given by 

dx'^ = dw"^ + dx^ + dy"^ + dz^ , (5) 

and the full four- dimensional space-time line element reads 

ds^ = c^dt^ - R^{t)dx'^ = ai.'nfldrf -dx^] . (6) 

Here R{t) = 0(77) denotes the cosmic scale factor as a function of cosmic time t 
respectively conformal time t] [dr] = cdt/R{t)). The metric © has the most general 
form assuming that the Universe is homogeneous and isotropic (cosmological principle). 

The condition w'^ + x'^+y'^ + z'^ = 1 can be used to eliminate one spatial coordinate, 
w say, and one is left with the three space variables {x,y,z). In the following, we shall 
also use the three angular variables (x, 0, 0) defined by 

w = cos X , a; = sin X sin 6' cos , 

y = sin X sin 6* sin , z = sin x cos ^ , (7) 

with 0<x^^5 0<^^<7r, O<0< 27r. It is seen from ((Zj) that (x, y, z) describe for 
a fixed value of x an iS^-sphere with radius sinx and volume V^^ {x) = ^(2x ~ sin2x), 
where 6 and play the familiar role of polar and azimuth angle. The 5^-sphere with 
radius sinx can be considered as a cross section of the unit 5^-sphere. Varying x 
between and vr, one obtains an infinite sequence of cross sections (=2-spheres) of S^, 
whose radius grows from to 1 (for < x ^ f )) and then shrinks again to zero (for 
f < X < ^); which provides a useful way to visualize the 3-sphere. 
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With (0) and ((7)), the hne element on can be rewritten as 

which has to be inserted into eq. ©• 

For a study of the large-scale anisotropics in the CMB produced by scalar 
perturbations, we need the (regular) eigenmodes of the covariant Laplacian A on S^, 
i. e. the solutions of the Helmholtz equation (A + £''^^)\E^'^^ = describing the vibrations 
on . Using the coordinates (x, ^, 0) and the metric (jH)) and separating the variables, 
the normalized solutions read iSl Hi ESI HEl El (^5 G N; / = 0, 1, . . . , /3 - 1; m = 
-/,-/ + !,...,/) 

m^{x) = Rpi{x)Yi^{eA) ■ (9) 
Here Yim{0,(f)) are real spherical harmonics on the unit sphere S"^, and the "radial 
functions" R/3i{x) are given by 

RM = Af,i {sin xYC^ft-Uicosx) , (10) 

where C^{x) are the Gegenbauer polynomials (which can also be expressed in terms of 
the associated Legendre functions P^{x)). The normalization factor in pUj) 



- ^ (/? + /)! ^^^^ 

follows from the orthonormality relations 

d^^{S)^flm{S)^?l'm'{^) = Sp,3' Sw 6mm' , (12) 

/ dx sin^ xR/siix) Rp'iix) = ^pp' ■ (13) 

The volume element on reads (with the standard expression dVL = sm6d6d(j) on S"^) 
dfi{x) = sin^xdxd^ ■ (14) 
Notice that the radial functions Rpiix) are even or odd (depending on whether j3 — l — l 
is an even or odd integer) about the equator x = f ? since they satisfy the periodic 
boundary condition 

Rpi{T^-x) = i-lf-'-'RM ■ (15) 

Since is compact, the vibrational modes are discrete, 

Ef = P'-l , (16) 

where /? G N plays the role of the wave number. The fundamental mode of has an 
eigenvalue of Ef^ = corresponding to the constant wave function "^ff^Q = l/W^, 
where V^^ = J^g d^{x ) = 27i'^ is the volume of the unit sphere S^. Since the eigenvalues 
of do not depend on / and m, the modes are degenerate, and the multiplicity of the 
mode j3 is given by 

r^\p) = J](2/ + l) = . (17) 
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If N'^ [k) denotes the number of the vibrational modes on with wave number smaller 
than or equal to A; G N, we obtain 

N^\k) = yr^'\l3) = -k' + -k^ + -k='^k^ + Oi^) 
^ ' 3 2 6 67r2 ^ ' 

13=1 

in agreement with Weyl's asymptotic law for k oo. 

We now turn to the multi-connected spherical spaces S^/T, T C S0(4), possessing 
the volume V{S^/T) = V^^ /N, where N denotes the order of the group T. (For 
reviews, see jlOl EU HHJ HHl ED])- To define the discrete fixed-point free subgroups 
r C SO (4) of the isometries of S^, it is advantageous to work with Hamilton quaternions 
q := w + xi + yi + zi], {w, x, y, z) G M^, which are spanned by the 4 basis quaternions 
{l,i,j,ij}, with the multiplication defined by i^ = = —1, ij = — ji plus the property 
that i and j commute with every real number. The multiplication of quaternions is 
associative, but not commutative. If q* := w — xi — y] — zi] denotes the quaternion 
conjugate to g, the norm of q is defined by |gp := qq* = w"^ + x"^ + y'^ + z"^ . Quaternions 
with norm 1 are called unit quaternions. For a unit quaternion g, there exists an inverse 
q~^ which is simply given by q~^ = q*. It is obvious that the unit 3-sphere can 
be identified with the multiplicative group of unit quaternions. The distance d{qi,q2) 
between any two points qi and q2 on is given by cos d{qi, ^2) = WiW2+XiX2+yiy2+ZiZ2 
and is a point-pair invariant, i.e. (^(751, 7^2) = d{qi,q2) V7 G S0(4). 

The group S0(4) is isomorphic to x 5^/{±(l, 1)}, the two factors corresponding 
to the left and right group actions. In the following, we are only interested in subgroups 
r G S0(4) which lead to homogeneous manifolds S^/T. Such groups possess only 
so-called right-handed Clifford translations. A quaternion 7 G F corresponds to a 
Clifford translation, if it translates all points qi,q2 G by the same distance 
i.e. (i(gi,7gi) = d{q2,'yq2) = X- The right-handed Clifford transformations 7^ G F 
in are defined by left-multiplication of an arbitrary unit quaternion g G 5^ by 
7fc := «fc + + Cfcj + i- e. 

q ^ lkq= (flfc + bki + c^j + 4ij) ■ {w + xi + yj + zij) 
= {akW - bkX - Cky - dkz) 
+ ipkW + akX - dky + Ckz) i 

+ {ckW + dkX + aky - bkz) j (18) 

+ {dkW - CkX + buy + akz) ij . 

The right-handed Clifford translation acts as a right-handed corkscrew fixed-point 
free rotation of . 

For our numerical computations of the eigenfunctions, it is convenient to represent 
the action ()18|) of a group element 7^ G F by the following orthogonal 4x4 matrix 





ak 


-bk 


-Ck 


-dk ^ 




bk 


ak 


-dk 


Ck 




Ck 


dk 


ak 


-bk 




dk 


-Ck 


bk 


ak ) 
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which acts from the left in an obvious way on the basis quaternions (1, i, j, ij). 
The following groups lead to homogeneous manifolds S^/T ^Ul US} I49j : 

• The cyclic groups Z„ of order n {n > 1). 

• The binary dihedral groups of order 4m (m > 2). 

• The binary tetrahedral group T* of order 24. 

• The binary octahedral group O* of order 48. 

• The binary icosahedral group /* of order 120. 

Let us consider the binary icosahedral group /* in more detail consisting of = 120 
group elements. They are generated by the two right-handed Clifford translations jSU] 

71 : g t-> 7i g := jg = -y + zi + wj - xij 

72 : g ^ 12 q 

with the unit quaternion 72 := f + + |j corresponding to the angles {X2,d2,(p2) = 
(|, |, arctano"), where a = is the golden ratio. The fundamental cell (Dirichlet 

domain) of the binary polyhedral space V := S^/I* is a regular dodecahedron and is 
known as the Poincare dodecahedral space, which is made of 12 pentagons, 30 edges and 
20 corners. 120 dodecahedra tessellate the unit 3-sphere S^. The shortest translation 
distance in V is given by d{q,'j2Q) = X2 = f , while 71 translates all points already to 
the much larger distance Xi = f • 

The vibrations on the general homogeneous spherical 3-spaces Ai^ = S^/T are 
determined again by the Helmholtz equation on S^, but now with the periodic boundary 
conditions corresponding to the defining fundamental cells (polyhedra). Thus the 
eigenfunctions "^^^{x) on can be expanded into the eigenf unctions \E''^^(a;): 

13-1 I 

^f'\^) = E E C,in.iM')^f,^{x) , (20) 

1=0 m=~l 

with E^^ = 13"^ - j3 e i = 1, . . . , r-^^(/3). It is important to notice that for 
a given manifold Ai"^ the wave numbers f3 do not take all values in N. The allowed 
j3 values together with their multiplicities r^^{P) are, however, explicitly known |51j . 
The full information of the non-trivial topology of A^'^ is contained in the real expansion 
coefficients {C/simi-^^)} which by virtue of (fT^ are given by 

CpimiM') = I d^,{x)^%^{x)^f\x) . (21) 
They satisfy the normalization condition 

/3-1 I 

E E i^hmiM')Y = N (22) 

1=0 m=-l 

due to the orthonormality relation 

/ d^^ix)^f'\x)^f'^'{x) 

= "^^^^^ "^f"^^^ ^^f'^'i^) = SpfS' ■ (23) 
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Summation of the relation (j22p over the degeneracy index i yields 

E E E = Nr^'m . (24) 

Z=0 m=—l 1=1 

In section El we shall determine the coefficients ^^;^(A^'^) numerically up to some wave 
number cut-off /^max- Here we state only the following Conjecture for homogeneous space 
forms (0 </</?- 1) 

m=—l i=l 

for which a numerical check will be discussed below. Let us also note that in our 
application to the CMB anisotropy, we only need the modes with /3 > 3 (if they exist), 
since the values /? = 1, 2 correspond to modes which are pure gauge terms 



3. Numerical determination of the eigenmodes 

For the homogeneous 3-manifolds Ai^ = S^/T, which we are considering here, the 
eigenvalues of the Helmholtz equation are given by = - 1, /3 e N, and even 

their multiplicities r^^{P) are explicitly known [HI]. As already mentioned in section |2l 
the wave number /3 does not assume all values in N, and thus the multiplicity r^^{P) 
refers only to the /3 values which appear in the spectrum. E. g. in the case of the Poincare 
dodecahedral space T>, there exist no eigenmodes for even /3-values, and, furthermore, 
also among the odd values jS G 2N + 1 there are finitely many gaps, i. e. (3 values for 
which no eigenmode exists. Explicitly, one has [HI] 

E^ = P'-l (26) 

with p e {1, 13, 21, 25, 31, 33, 37, 41, 43, 45, 49, 51, 53, 55, 57} U{2n + l,n> 30} and 

The fact that the eigenvalues respectively the wave numbers (3 for a given 
manifold are explicitly known, facilitates greatly the numerical computation of the 
eigenfunctions '^{S), since it saves us from the time-consuming numerical search 
for the allowed /^-values as it has to be carried out, e.g. for hyperbolic manifolds in 
the case of a negatively curved universe I54j . Furthermore, since the expansion ()2()j] 
in terms of the eigenfunctions involves for a given wave number only a finite basis 
of dimension (see eq. (|T7|l). the number of the real coefficients {Cpimi-^^)} is also 
restricted to (3"^ for a given eigenfunction with wave number f3 and degeneracy index i. 

The eigenfunctions have to satisfy the fundamental periodicity condition 

^f''{lkq) = ^f^'{q) , yqeS' , V 7. G F , (28) 

where 7^ denotes a group element of the group F C SO (4) which defines the manifold 
M'^ = S^/T. We will use the condition ^ in a collocation algorithm. Thus, the 



r^(/5) = P 



P-1 
10 



+ 



P-1 



+ 



P-1 



Poincare Dodecahedron 



10 



condition (j2H)) will be imposed on L > randomly distributed points q on S^. 
To generate such random points, we use instead of the coordinates {w,x,y,z) the 
coordinates (f , a, 7) defined by 



w = y/1 — 2v cos a 

X = cos 7 ^^gN 

y = \/2v sin 7 



z 



\f\^^^ sin a 



with V G [0, |] and a, 7 G [0, 27r]. The line element (0) takes then the form 

dx^ = TTT-. ^ + (1 - 2f ) rfa^ + 2vd-i'^ , 

2f (1 — 2v) 

which gives the volume element d^ = dvdad'j. The random points q on are now 
generated by randomly choosing the coordinates (f,a,7) and then calculating from 
pnjl the coordinates {w,x,y,z) of the corresponding unit quaternion q. In this way, we 
generate L random points {gi, g2, • • • , ^l}- An arbitrary group element 7fc G F maps the 
point qr to the point := Ikqr on S^, where in the following computation the group 
elements 7^ run over r\{l}. The action of 7^ will be represented by the matrix M^, 
see eqs. (|THtfTTI|) . After having generated in this manner L mirror points {qi, (?2, • • • , ^l}, 
one imposes the condition at all those L points. Inserting the expansion (j^ . one 
finally obtains the following over-determined system of equations (r = 1, . . . ,L > (3'^) 
for the coefficients {^pimi.-^^)} ■ 

/3-1 / 

E E [Km - Kln^i^r)] e^lmiM') = . (30) 
1=0 m=-l 

The system (|Hn|l is for a given /5-value of the form J2s=i ^rs^s = 0, r = 1, . . . , L, where 
the (L X /3^)-matrix B = B{[3) is given by the square bracket in eq. (jHUj) and is thus 
completely determined in terms of the eigenfunctions on the covering space . The 
vector ^ = ^{(3) has the components = ^pimi-M^), i-e. s denotes the multi-index 
{00, 1 — 1, 10, 11, . . . , /3 — 1 /5 — 1}. The system flSUj) is numerically solved by using the 
singular value decomposition. One then obtains (P) linearly independent solutions 
for the expansion coefficients {Cpimi-^^)} which are distinguished by the index i. In a 
final step, the coefficients are normalized according to the condition (|22|l . 

Using this method, we have computed the expansion coefficients for F = /*, i. e. 
for the Poincare dodecahedral space "D = S^/I*, for j3 < 155 comprising the first 10 521 
eigenfunctions. The normalized coefficients have been used to check the conjecture ()25j) 
which has been found to hold within a numerical accuracy of 13 digits. In addition, 
we have computed the coefficients {Cpimi-^^)} checked the conjecture (^Hjl for 
/3 < 33 for some cyclic groups, some dihedral groups, the binary tetrahedral and binary 
octahedral group and have found again excellent numerical agreement. However, for 
the inhomogeneous lens spaces [50] L(12,5) and L(72,17), we have found that the 
conjecture does not hold. We thus conclude that the conjecture is only valid 
for the homogeneous 3-spaces. 
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4. The angular power spectrum ST^ and the correlation function C(^9) 

In order to determine the angular power spectrum 6T^ and the temperature correlation 
function C{i!}), one has to compute the temperature fluctuation ^{n). This is done 
using the Sachs- Wolfe formula which is given within the tight-coupling approximation 

by (c = 1) 



^(n) = / d^k 



+ 2 d^k / rfr^^iiiZvI/^,(r(77),^,0) , (31) 

where the fc-integration has to be replaced by a summation over j3 in the case of a discrete 
spectrum. (Here ?7sls denotes the conformal time at recombination corresponding to a 
redshift 2;sls = 1089.) are the eigenfunctions on Ai^, and $g(?7) are the expansion 
coefficients of the metric perturbation $, see eg. (p?^ below. {r{ri),9,(f)) denote the 
spherical coordinates of the photon path in the direction n, where we now assume that 
the observer is at the origin of the coordinate system, i. e. at (xobs, 6'obs, 0obs) = (0, 0, 0). 
Here, %^(?7) is the expansion coefficient of the relative perturbation in the radiation 
component, and Vj:^{ri) is the expansion coefficient of the spatial covariant divergence 
of the velocity field of the tightly coupled radiation-baryon components. The method by 
which the quantities occurring in (jHT|) are numerically computed, is described in detail 
in Section 2 of [SEj- 

The ordinary Sachs- Wolfe (SW) contribution to the temperature fiuctuation 6T{n) 
is a combination of the gravitational redshift due to the gravitational potential 
^{ri,T,9,(j)) at the surface of last scattering (SLS) and the intrinsic temperature 
fiuctuation due to the imposed entropic initial conditions. Since the SW 
contribution, i.e. the first two terms in the brackets in (see also eq. (jHSI) below), 
gives a dominant contribution to the CMB anisotropies at large scales, we first present 
an exact analytic expression for the expansion coefficients {aim}, the multipole moments 
{Ci} and the temperature correlation function C{'d) for the Poincare dodecahedral 
universe taking only the scalar perturbations into account due to the SW effect. 

For an energy-momentum tensor T^j, with Tij = for i ^ j and i,j = 1,2,3, the 
line element can in conformal Newtonian gauge be written as [22] 

ds"^ = a^{r]) [(1 + 2$)dr/2 - (1 -2$)|rfxp] . (32) 

Denoting hj h = n{6, (p) the unit vector in the direction in which the photons are 
observed, the SW formula reads 

7^ = HVSLS, T-SLS, 0, 0) + -5^(r/sLS, T-SLS, 0, 0) (33) 

- ^*(%LS,1-SLS,^,0) (34) 

with TsLs '■= Vo ~VsLS- (Here ?7o denotes the conformal time at the present epoch.) The 
SW contribution (jH^ can be approximated as in (jH^ for modes well outside the horizon 
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at the time ?7sls of l^^st scattering, i.e. for modes with P -C The approximation 

is used in the following derivation of the analytic expressions for Cf^ , eq. (jS}, and 
C^^{^), eq. (jlH). In the actual numerical computation, we use (pTH) as discussed above. 
For the dodecahedron V, the metric perturbation can be written as an expansion in the 
eigenfunctions of V 

Hv,r,e,<t>) = Yl' E '^M^)^?'(^'^''^) • (35) 

/3>13 i=l 

Here the prime in the summation over the modes /3 indicates that f3 G 2N + 1 does 
not assume all values. The functions $^(77) determine the time-evolution and will be 
factorized 

$^(ry) = ^i{0)gpiv) (36) 

with 5'/3(0) = 1. The functions gpir]) do not depend on i, since the associated differential 
equation depends only on the eigenvalue which is independent of i. The initial values 
$^(0) play the role of the primordial fluctuation amplitudes and are Gaussian random 
variables with zero expectation value and covariance 

$^(0)$J,(0)) = 6fs^> 6,, P^{(3) . (37) 

Here Pq>{l3) denotes the primordial power spectrum which determines the weight by 
which the primordial modes (3 are excited, on average. The average (. . .) in ()37|) denotes 
an ensemble average over the primordial perturbations. In the following, we shall assume 
the scale-invariant Harrison-Zel'dovich spectrum 



The only free parameter in (jH8|) is the normalization constant a which will be determined 
from the CMB data. 

Inserting the expansion into the Sachs- Wolfe formula (j34|) and using the 
eigenfunctions of the dodecahedron (see eqs. (j^Uj) . and ((H)), we obtain for the 
expansion coefficients of the SW contribution (/ > 2) 

af:^m = 3E' E'^M^SLs)a/™(I^)i?/3/(rsLs) , (39) 

f3>l 1=1 

which in turn yields with the definition Q and the Gaussian hypothesis ()37|1 the mean 
values of the multipole moments (Z > 2) 



9 

I3>1 



I r-^ifi) 



2/ + 1 . , 

m=—l 2=1 



<(^SLS) . (40) 



At this point we can use our Conjecture (j25|l with = 120 and arrive at the final form 
CfW(D)^fE'^^^^4(*.s)4,(--) . (41) 

/3>/ ^ 
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The last expression shows that the lowest multipoles are suppressed due to the discrete 
spectrum of the vibrational modes of T>. Since the lowest modes are at (3 = 13, 21, 25 
and 31, the mode sum starts at /5 = 13 for the multipoles / = 2, . . . , 12, while for the 
multipoles / = 13, . . . , 20 it starts only at f3 = 21, for the multipoles / = 21, . . . , 24 only 
at P = 25, and for the multipoles / = 25, . . . , 30 only at /3 = 31. 

From (PT|) we can compute the SW contribution to the correlation function ^ for 
the Poincare dodecahedron V 

C--(tf) = H^'!fM^,J(,,,,)C'(ftrf) , (42) 

/3>13 ^ 

1 /2 

where the angular dependence is completely given 

13-1 

C{P,{}) := ^(2/ + l)4^(rsLs)PKcos^) 

1=2 

13-1 

= 5^(2/ + l)4,(sinrsLs)'' [^^.^^.(cosrsLs)]' ^^//^(cos^) (43) 

1=2 

after having inserted the radial function (|10|). With the help of the addition theorem 
for the Gegenbauer polynomials (see p. 223 in ^^), one obtains for the mean value of 
the SW contribution to the correlation function for the Poincare dodecahedron universe 



20 r^i(3)P^iP) 2, , sin(/37(t9)) 



3vr2^^3 P ^^^'""'^ sin7(^) 

_ ^c<sw _ A^swcos^? , (44) 

where the relation C^_^(cos7) = ^^^-^ has been used. The angle 7 = 7('t9) G [0, 2rsLs] 
is given by 

7(7?) = arccos ( cos^ tsls + sin^ tsls cos -(9) , (45) 
and the monopole and the dipole moment read 

/3>13 ^ 



sm Tsls 



(46) 



Cr = ^sm2(rsLs)2^ ^2/^2 ^g(^SLs) [C2_2(cosrsLs)] • (47) 

/3>13 ^ ^ 

In figure El we show as a dashed curve the evaluation of eq. (jlH) for f2tot = 1.020. 
The solid curve shows the full result including the integrated Sachs- Wolfe and Doppler 
contribution. The grey band corresponds to the Icr deviation computed from 500 
simulations. Although the curve corresponding to the ordinary Sachs- Wolfe contribution 
lies within the la band of the full result, the other two contributions cannot be neglected 
as discussed in more detail below. 

In our computations of the multipole moments Ci for a given primordial realization, 
we use all eigenmodes with f3 < 155. Above this cut-off, we use in the mode summation 
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Figure 2. The SW contribution l|44|l of the temperature correlation function C(i}) is 
shown as a dashed curve for the dodecahedral space V for fitot = 1-020. The solid curve 
shows the full result including the integrated Sachs- Wolfe and Doppler contribution. 
The la deviation is shown as a grey band. 

up to (3 = 1501 the conjecture ()25|) which describes the mean asymptotic behaviour. If 
only the modes below (3 = 155 are used, the multipoles can only be calculated up to 
/ ^ 15. It is thus important for the calculation of (^(-i?) that the conjecture allows us to 
include higher multipoles. 

Let us now come to the discussion of the CMB anisotropy in the Poincare 
dodecahedral space, where, for simplicity, we assume that the dark energy is given 
by a cosmological constant with Q\ := Q^ot — ^rad — ^mat- We have calculated the CMB 
anisotropy for different values of the cosmological parameters. Here we present only the 
results for fibar = 0.046, fi^at = 0.28 and h = 0.70. 

20000 

s 

15000 
10000 
5000 


1.01 1.015 1.02 1.025 Otot 1.03 

Figure 3. The values of the S{p) statistic (|48ll for the median values of C{i3) for the 
dodecahedral space for p — 60° (full curve) and p — 20° (dotted curve) in dependence 
on Qtot {h = 0.70, f^mat = 0.28 and f^bar = 0.046). 

To quantify the observed surprisingly low CMB anisotropy 4 large scales, the 
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S statistic 

/cos p 
\C{^)\^d cos^ (48) 

is discussed for the first year WMAP data in [2H] for P = 60°, and it is found tliat only 

0. 3% of the simulations based on the concordance model ri have lower values of 5'(60°) 
than the observed value S'(60°) = 1644. Somewhat higher values of S'(60°) are obtained 
using other statistical methods and other sky masks in [HSl, but they are nevertheless 
surprisingly low. In figureElwe show the values of S'(60°) (full curve) and 5(20°) (dotted 
curve) using the median values of C for the dodecahedral space V in dependence on 
f2tot- One observes that the models with i^tot = 1-017 . . . 1.020 give the lowest values for 
the large scale anisotropy. 

Here and in the following calculations, the overall amplitude of the CMB anisotropy, 

1. e. the normalization constant a in fl38|) . is fitted to the Ci values of the first year WMAP 
data in the range I e [20, 45]. (Thus, the values cannot be directly compared to the ones 
in [33] , where the amplitudes are scaled such that they match the C4 value of WMAP 
exactly.) 




30 60 90 120 150 I? 180 30 60 90 120 150 1!) 180 




Figure 4. The temperature correlation function C(i?) for the dodecahedral topology 
for four values of f^tot for h = 0.70, f^mat = 0.28 and f^bar = 0.046, the WMAP first 
year observation and the pl-concordance model (see description in the text). 



In figureElthe temperature correlation function C{'d) is shown for four values of fitot- 
The WMAP curve obtained from the LAMBDA home page ^http: //lambda. gsfc.nasa.gov 
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is shown as a full curve. The mean value of C{-&) for the pl-concordance model from 
the same home page is displayed as dotted curve. The la deviation computed from 
3000 HEALPix simulations is shown as a shaded region. Finally, the mean value of 
C for the dodecahedral topology is shown as a dashed curved and the corresponding 
la deviation is the grey region obtained from 500 simulations. In figure E^), the value 
preferred in [23], i.e. fltot = 1-013, is presented. One observes a general low anisotropy 
on large scales, however, a more than la deviation from the WMAP values occurs in 
the range G [60°, 90°]. The WMAP data are much better described by models with 
^tot = 1-018 and Qtot = 1-019 (figures IHd) and Eh)) as is also revealed by the better 
values for the S statistic (compare figure E}- Only on the largest scales ^ > 170°, the 
models in figures EJd) and Et) do not reproduce the large negative correlation obtained 
from the WMAP data. But without this blemish, the dodecahedral space describes in 
the considered Qtot range the observations down to the smallest scales, since the WMAP 
curve lies completely within the grey la region. Finally, figure Eli) shows a model with 
^tot = 1-026, where the S statistic already deteriorates and consequently the agreement 
with the WMAP curve is also worse. 




Figure 5. The angular power spectrum ST^ for the dodecahedral topology (circles) 
for four values of fltot for h = 0.70, f^mat = 0.28 and f^bar = 0-046. The WMAP data 
are shown as diamonds with error bars, and the pl-concordance model is represented 
by the full curve. 



In figures EK) to Eli), the angular power spectrum 6T^ = l{l + 1)Ci/(2ti) is shown 
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for the same values of Qtot as in figure EJ The WMAP data are represented by diamonds 
together with the la error bars which do not include the cosmic variance. The mean 
values of ST^ of the pl-concordance model (full curve) do not show the suppression of 
power at low values of I, but in contrast, increase there. The mean values of 6T^ for the 
dodecahedral space exhibit the suppression of power as observed by WMAP. A good 
agreement is observed exactly for those values of Qtot favored by the S statistic and 
the correlation function C{i}). This behaviour is also revealed by figure IHl which shows 
the dependence on Qtot of the mean values of 6T2, ST^ and (5T| for the dodecahedral 
space V. The lowest values of the quadrupole 6T2 are obtained near f^tot = 1-017. 
Furthermore, at fitot = 1.017 the next multipole 6T^ matches the value obtained by 
WMAP almost perfectly. 




' ' 

1.01 1.015 1.02 1.025 Otot 1.03 

Figure 6. The mean values of 6T2 (solid curve), 5T^ (dashed curve) and ST^ (dotted 
curve) are shown for the dodecahedral space I? as a function of fitot- The lower 
horizontal line indicates the (5r| value of WMAP and the upper one the corresponding 
value of 6Tl 

The probability distribution of ST2 and ST^ obtained from sky simulations is for 
the pl-concordance model much broader than for the dodecahedral space. To emphasize 
this fact, figure d displays the distribution of 6T2 and ST^ for f2tot = 1.018. The 
corresponding WMAP values are indicated by the vertical lines together with the la 
error (vertical dark grey band). The shaded histogram shows the wide distribution 
of 3000 HEALPix simulations for the pl-concordance model. In contrast, a much more 
pronounced peak is revealed by the corresponding distribution of the dodecahedral space 
(light grey histogram). Furthermore, the peaks for 6T2 and 6T^ lie both within the la 
error of the WMAP observation. 

5. The circles-in-the-sky signature 

A multi-connected space can reveal itself by the so-called circles-in-the-sky signature 
Consider an observer located at a given point x who measures the fiuctuations of 
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Figure 7. The probability distribution of ST^ and ST^ in [^K^] for fltot = 1.018, 
h — 0.70, r^mat = 0.28 and f^bar = 0.046 for the dodecahedral space (hght grey 
histogram) . The wide distribution of the pl-concordance model is shown as the shaded 
histogram. The WMAP values together with their la error are indicated by the vertical 
line and the vertical grey band, respectively. 



the CMB. The contribution to the CMB fluctuations due to the ordinary Sachs- Wolfe 
(SW) effect arises on a sphere around this observer having the radius tsls, i-e. at the 
so-called surface of last scattering (SLS) w. r.t. this observer. In the universal cover 
there are copies of this observer which have their own SLS's which are centered at the 
copies of the observer. Those copies whose distance to the observer at x are smaller 
than 2rsLs possess SLS's which intersect the SLS of the observer sitting at x. Since the 
intersection of two spheres is a circle, both the observer and its copy observe along this 
circle the same temperature fluctuations due to the ordinary Sachs-Wolfe effect. Since 
the observer and its copy have to be identified, there are two matching circles having 
the same radius on the sky and having an identical ordinary Sachs- Wolfe contribution. 
The integrated Sachs- Wolfe (ISW) contribution arises on the photon path from the SLS 
to the observer, see eq. (j!?T|l . Since the two paths are not identified, their contribution is 
different on the matched circles. The two Doppler contributions are also different, since 
their velocity projections are different. 

In the case of a flat toroidal universe, these different contributions to the 
temperature anisotropy are investigated in [SHj with respect to the circles-in-the-sky 
signature, and it is found that the Doppler and the integrated Sachs- Wolfe contributions 
degrade the signature. That these two contributions play an important role also in the 
case of the dodecahedral topology, is demonstrated in figure |H1 It shows the three 
competing contributions as computed from eq. ()31|) for the simply connected (figure 
ISt)) and for the dodecahedral space V (figure |Hb))- Here we use the conjecture 
up to /3 = 4001. For the SW contribution, see eq. ()41|) . The most obvious difference 
between V and iS^ is the suppression of power for small multipoles 6Ti in the case of 
V in comparison with S^. In the case of the dodecahedral space V, the loss of power 
is so severe that for the quadrupole, the ISW dominates over the ordinary Sachs- Wolfe 
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contribution. A further distinction between and V is that STi oscillates as a function 
of / in the case of the multi-connected space even for the mean values of 6Ti which 
are considered here. With the exception of the smallest multipoles, the ordinary Sachs- 
Wolfe contribution dominates up to /< 30. However, even in this range the Doppler and 
the ISW contribution provide a significant fraction of the total temperature fluctuation. 
The ISW is most important for I < 10, whereas the Doppler contribution increases for 
/ > 10 to even larger values than the ISW, such that both contributions give together a 
large perturbation to the pure circles-in-the-sky signature. 




To emphasize this point, figure M shows the temperature fiuctuation along two 
matched circles parameterized by the angle defined in fl49j) below. One observes 
in figure ^) that the total temperature fiuctuation 6T is not matched perfectly well. 
However, some rough similarities between both 6T curves are nevertheless visible. 
Whether such similarities occur sufficiently frequently by chance in the simply connected 
model, such that the circles-in-the-sky signature is swamped, is discussed below. The 
different contributions to 6T are separately shown in figures Eb), Efc) and Eli)- Even the 
ordinary Sachs- Wolfe contribution (figure inb)) does not provide a perfect match which 
is due to the fact that the dipole contribution has been subtracted from the sky map 
analogously as it is done on the observational side. This destroys the perfect agreement 
due to the SW contribution. Conversely, if one would know the topological structure of 
the Universe, this would offer an opportunity to determine the dipole contribution due 
to the primordial fluctuations which is usually superseded by the Doppler shift due to 
our local motion. The Doppler and the ISW contributions are completely different for 
the two circles as seen in figures Ot) and Oil), respectively, for reasons described above. 

As just discussed, for cosmological models near fitot — 1-02, the ordinary Sachs- 
Wolfe contribution dominates for / < 30. The value / = 30 corresponds to a scale 
9 ~ — 6°- Thus the circles are blurred on scales below ~ 6°. On larger scales, the 
integrated Sachs- Wolfe and the Doppler contribution lead to a modulation of matched 
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Figure 9. The temperature fluctuation ST is shown along two matched circles for 
Qtot = 1-019, h = 0.70, 17mat = 0.28 and f^bar = 0.046. Panel a) displays ST of the 
total temperature fluctuation along the two circles, obtained from the ordinary Sachs- 
Wolfe (SW), Doppler (Dop) and integrated Sachs- Wolfe (ISW) contribution, which are 
shown separately in panels b), c) and d), respectively. 



circle structure. 

For a quantitative search of matched circles in microwave sky maps, the quantity 

is introduced in (In IBj, E is called S, a name we avoid in the following in order 
to not confuse the reader with the previously discussed S statistic.) Here STa{(p) and 
STb{(f)) are the temperature fluctuations along two circles on the SLS with the same 
radius, and () := J^^ d(f). The angle p describes the relative phase between the 
temperature fluctuations on the two circles and the plus/minus sign in the nominator 
allows for orientable as well as for non-orientable manifolds. The S variable takes 
the value E = -|-1 in the case of a perfect correlation and S = — 1 for anticorrelated 
temperature fluctuations. 

In the E test is applied to the first year data from WMAP and it is found 
that Emax = niaXp{E(p)} does not show any spikes which would reveal a matched circle 
pair. The search is carried out for nearly back to back circles having a separation larger 
than 170° and a radius larger than 25°. For 1.0091 < Qtot < 1-0252 (for h = 0.70 
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Figure 10. The radius of the six matehed circles is plotted as a function of fitot- 

and l^mat = 0.28) the Poincare dodecahedron has in total 6 pairs of circles which are 
back to back, i. e. are separated by 180°. The radius of these circles increases with 
increasing f^tot as shown in figure ITUl Above f^tot = 1-0252 there are at first 16 pairs, a 
number which increases further with increasing Qtot, and below Qtot = 1-0091 there are 
none. The conclusion in [SHI is thus that the topology of the Poincare dodecahedron is 
excluded. However, in jHO] a hint towards six pairs of circles having a radius of 11 ± 1° 
is found which would correspond to i^tot = 1-0094 for our choice of h and fimat- 

Now we would like to discuss whether the Doppler and ISW contributions are 
significantly large enough in order to contaminate the SW contribution such that the 
circles-in-the-sky signature is hidden. We generate 500 sky maps for the dodecahedral 
space V using all modes up to /3 < 155 and compute for each of the six pairs of matched 
circles the value of S^ax- In this way we get 3000 values for Smax belonging to circle pairs 
including the SW, Doppler and ISW contribution. In addition, we obtain 3000 values 
for where p is not varied but instead is fixed to the value p = 36° determined by 
the topology of V. Due to the Doppler and the ISW contribution, the values of Smax are 
usually larger than the correct values S(p) for the matched pair of circles having p = 36°. 
This emphasizes that the contamination by the Doppler and the ISW contribution is so 
large that the Smax statistic does in general not find the correct value of p. Furthermore, 
we compute also the corresponding values of Smax for 500 simulations for the simply 
connected universe with the same cut-off in (3 and with the ring positions found in 
the dodecahedral space V. Of course, the Emax values of do not correspond to a 
matched circle pair since has none. 

In figure [TTl we compare the probability distribution of S^ax for and V for 
the two values Qtot = 1.011 and Qtot = 1.020. These values are chosen in order to 
demonstrate the effect due to the size of the radii of the circles which are 25.4° and 
49.8° for Qtot = 1-011 and Qtot = 1-020, respectively. The probability distributions 
shift to smaller values of S^ax with increasing radii. In figure ^1 we again show the 
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probability distribution of Smax for as in figure [TTJ but display for V the probability 
distribution of S(p) with p = 36°. It is seen that the distribution for V has been 
shifted to smaller values of S, and thus the overlap of the distributions for V and 
has increased. 
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Figure 11. The probability distribution of Emax for (shaded histogram) and 
V (grey histogram) for iltot ~ 1.011 in panel a) having a radius of 25.4° and for 
^^tot — 1.020 in panel b) having a radius of 49.8°. The other cosmological parameters 
are h = 0.70, ilmat = 0.28 and flhar = 0.046 as usual. 
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Figure 12. The probability distribution of Smax for (shaded histogram) as in 
figure [TTl and of ^(p) for V with p — 36° as determined by the topology of V (grey 
histogram) for fJtot = 1.011 in panel a) having a radius of 25.4° and for fitot = 1.020 in 
panel b) having a radius of 49.8°. The other cosmological parameters are as in figure 

im 

Analysing a sky map which possesses no pairs of circles, one would obtain for 
randomly chosen circle pairs a distribution corresponding to the shaded histograms in 
figures ITT] and IT^ Now the aim is to find the six pairs of circles due to the dodecahedral 
topology. If there would be only the pure ordinary SW contribution, a value of Smax = 1 
would clearly stand out of the distribution due to S^. However, the Doppler and ISW 
contributions degrade this signal to the grey histograms in figures ^2 and ^| Thus 
the six pairs would probably yield values of Smax between 0.4 and 0.8 for f2tot = 1.011 
and between 0.4 and 0.7 for fltot = 1.020, and values of S(p) between 0.3 and 0.75 
for both values of f2tot, respectively. In the case of an analysis of a real sky map, one 
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would obtain many thousands of S^ax values which would not correspond to matched 
circles, i. e. having a probability distribution looking like that of S^. Embedded in this 
distribution, one would for the dodecahedron obtain six E^ax values corresponding to 
real matched circles. Because of the large overlap of both distributions, the signal of the 
six pairs of circles could easily be swamped. Thus there is a non-vanishing probability 
that the six pairs of matched circles can be overlooked in the sky maps. The search for 
matched circles is in addition made more difficult by the finite thickness of the SLS and 
the Sunyaev-Zeldovich effect, a possible systematic in the removal of foregrounds and 
further secondary effects which are not included in our analysis. 

6. Discussion and Summary 

Taking the WMAP value Qtot = 1-02 ± 0.02 as its face value, we conclude that the 
Universe is slightly positively curved and therefore necessarily spatially finite, since all 
spherical space forms possess this property. At present, the only way to decide on the 
particular topology realized in our Universe consists in studying a concrete space form 
and then compare its predictions with the observational data. 

In this paper, we produce evidence in support of the hypothesis that the spatial 
structure of the Universe could be given by the Poincare dodecahedral space V, provided 
that the total energy density parameter f^tot is in the range 1.016 ... 1.020. For the 
detailed comparison with the WMAP data, it was crucial that we could base our 
predictions for the CMB anisotropy upon a large number of vibrational modes of 
the Poincare dodecahedron. Since the eigenfunctions of the Poincare dodecahedron 
are not known analytically, they have to be computed numerically. Using the very 
efficient algorithm described in section 13 we could take in our calculations the first 
10 521 eigenfunctions into account. In addition, we put forward the Conjecture ()25j] 
for the homogeneous spherical space forms. Combining the numerically computed 
eigenfunctions with the Conjecture, we were able to calculate the CMB multipoles in 
a broad range, see figures El and |H1 covering the region of the large scale fluctuations. 
Fixing the cosmological parameters as fibar = 0.046 and flmat = 0.28, h = 70, in 
agreement with the WMAP data, there are only two parameters left: the amplitude 
a of the scale-invariant Harrison-Zeld'ovich spectrum and the density parameter 
fitot > 1 or, equivalently, the curvature radius Rq = -§^{i^tot — 1)^^^^. (The dark energy 
is identified with a cosmological constant with magnitude Q\ := Qtot — ^rad — ^mat-) 
The parameter a was obtained from a fit to the Ci values as determined by WMAP in 
the range I G [20,45]. (This differs from the normalization in [33], where a has been 
set to match the WMAP multipole with / = 4 exactly, since higher multipoles could 
not be calculated in [SHI due to the very limited number of available eigenfunctions.) 
Having fixed in this way all parameters except fltot, we studied the CMB multipoles and 
the temperature correlation function C{'&) for the dodecahedron as a function of f2tot- 
One should keep in mind that we restricted the parameter space to the above values. 
Thus we cannot exclude that other choices could lead to different values for Qtot due to 
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geometrical degeneracies. 

As a quantitative measure of the quality of our theoretical predictions, we calculated 
the S statistic (jlKjl for the angles p = 60° and p = 20°. The result is given as a function 
of Qtot in figure El and shows for this statistic the remarkable feature of pronounced 
minima at almost the same fltot values, Qtot G [1.016,1.020]. In particular, figure El 
shows that the Poincare dodecahedral universe leads for fltot in the above range to a 
suppression of the CMB anisotropy at large scales in agreement with the WMAP data. 
This is also demonstrated in figureEh), which shows the temperature correlation function 
C('i9), and in figures Efc) and IHb) showing the CMB power spectrum {fltot = 1-019 is 
used in these figures). Figure Et) demonstrates clearly that the dodecahedron model 
describes the WMAP data much better than the concordance model. At this point 
it is important to notice that the suppression of the CMB anisotropy at large scales 
occurring in the dodecahedron is neither obtained by fitting the theoretical curves to 
one of the first (suppressed) multipoles nor to the correlation function C . Rather the 
suppression is due to the finiteness of the Universe in combination with a small volume 
of the dodecahedron, ¥{1)) = ^ = 0.16449. . ., which is nicely illustrated in figure |H1 
where we compare the power spectra of the dodecahedron V and the simply connected 
iS^, respectively. The structures seen in figures El and |Hb) in the CMB power spectrum 
at multipoles / < 20, are mainly caused by the ordinary Sachs- Wolfe (SW) contribution 
(dotted curve in Eb)) whose analytic expression is given in eq. PT|) . The observed 
structures are a direct consequence of the discrete mode spectrum of V possessing large 
gaps at small wave numbers (3, see eqs. fEHj) and (j77|) . 

The particularly striking suppression of the quadrupole and octopole is also 
exemplified by figure which displays the probability distribution P{6T^) (Z = 2, 3) 
obtained from sky simulations. The probability distributions for the dodecahedron 
show pronounced peaks which overlap nicely with the la band of the WMAP data. In 
contrast, the concordance model predicts much broader distributions producing large 
values for 6T2 and (5T|, seen also in figure El 

Summarising the results discussed so far, it is fair to say that the Poincare 
dodecahedral universe provides a natural explanation for the strange suppression of 
power at large scales as observed by COBE and WMAP, if f^tot is in the range 
1.016 ... 1.020, which solves "the mystery of missing fiuctuations" . At this point we 
would like to emphasise that the good agreement of the Poincare dodecahedral universe 
with the WMAP data does by no means exclude that there are other manifolds which 
lead to similar conclusions. Rather, it may well be that there exists a "topological 
degeneracy" in addition to the well-known "geometrical degeneracy" (see e.g. [2^ I25j). 
Since the value f2tot = 1-02 ± 0.02 reported by the WMAP team depends on several 
priors and includes the la deviation uncertainty only, there still exists the possibility 
of the Universe having negative curvature. In fact, e. g. the hyperbolic Picard space 
[211 EZ| leads also to a suppression of the CMB fiuctuations at large scales in agreement 
with the WMAP data. 

There remains to detect a more direct signature of the dodecahedral topology. 
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For this purpose, we have studied the circles-in-the-sky signature in section 1^1 The 
dodecahedral universe possesses exactly 6 matched circles with radii 42.7° and 49.8° 
for f2tot = 1.016 and f2tot = 1-020, respectively. If there would be only the ordinary 
Sachs- Wolfe contribution, one would predict an almost perfect matching of the circles. 
However, as it is seen in figure |Hl there are important contributions from the integrated 
Sachs- Wolfe and the Doppler effect. The three contributions together destroy a perfect 
matching as is illustrated for a given circle pair in figure ^^). While the ordinary 
Sachs- Wolfe contribution displayed in figure Eb) shows the expected matching, the 
Doppler and the integrated Sachs- Wolfe contributions shown in figures Ot) and M^), 
respectively, differ appreciably over the full angular range, and thus lead to the poor 
matching demonstrated in figure IHt')- 

The quantity eq. P9|l . has been introduced for a search of matched circles 

in the microwave sky. Figures [TT] and [T21 display the probability distributions of Smax 
and S(36°), respectively. One observes broad distributions for both the dodecahedron 
V and the covering space S^. While the distribution for T> is peaked at higher values 
of Sinax than the distribution for S^, the two distributions have a large overlap. Due to 
the contamination by the Doppler and integrated Sachs- Wolfe effect, the distribution 
for V has no peak at Smax — 1, but rather gives for the most probable Smax values at 
which the 6 matching circles are expected to be observed the interval 0.4 . . . 0.8. 

In view of these complications, it seems that a search for matching circles in the 
first-year WMAP data is quite difficult, and thus the negative result reported in jH] 
does not yet allow to conclude that the dodecahedral model is excluded. We hope that 
the second-year WMAP data and the CMB measurements from the Planck satellite 
will improve the situation considerably such that we shall be able to decide about the 
topology of the Universe in the near future. 

Let us close with a quotation from Einstein, who wrote in 1919 in describing the 
Universe of beings in a "finite and yet unbounded Universe" [HD 102] : "The great charm 
resulting from this consideration lies in recognition of the fact that the universe of these 
beings is finite and yet has no limits." 
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